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ABSTRACT 

We report the discovery of very high energy gamma-ray emission from the direction of the SNR 
G54. 1+0.3 using the VERITAS ground-based gamma-ray observatory. The TeV signal has an overall 
significance of 6.8(T and appears point-like given the resolution of the instrument. The integral flux 
above 1 TeV is 2.5% of the Crab Nebula flux and significant emission is measured between 250 GeV 
and 4 TeV, well described by a power-law energy spectrum dN/dE ~ E"'" with a photon index F ~ 
2.39 ± G.23stat i 0.30sys- We find no evidence of time variability among observations spanning almost 
two years. Based on the location, the morphology, the measured spectrum, the lack of variability and 
a comparison with similar systems previously detected in the TeV band, the most likely counterpart 
of this new VHE gamma-ray source is the PWN in the SNR G54. 1-1-0.3. The measured X-ray to 
VHE gamma-ray luminosity ratio is the lowest among all the nebulae supposedly driven by young 
rotation-powered pulsars, which could indicate a particle-dominated PWN. 
Subject headings: gamma rays: general — ISM: supernova remnants — ISM: pulsar wind nebulae — 

pulsars: individual (J1930-f 1852, J1928+1746) — supernova remnants: individual 

(G54.1+0.3 = VER J1930+188) 
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1. INTRODUCTION 

In recent years, new insights into the most extreme en- 
vironments in the Galaxy have been obtained through 
the ve ry high energy (VHE, >100 GeV) gamma-ray 
band (Aharonian et al. 2005[ ). The majority of the 
>50 Galactic I'eV sources tound thus far are associated 
with supernova remnants (SNRs) or pulsar wind nebu- 
lae (PWNeFM The detection of these stellar remnants 
in VHE gamma rays constitutes a probe of non-thermal 
astrophysical processes. The acceleration of electrons or 
nuclei, above 1 TeV, in these sources leads to the genera- 
tion of gamma rays. Relativistic electrons can produce 
gamma rays by non-thermal bremsstrahlung, and by 
inverse-Compton (IC) scattering on ambient microwave, 
IR, or optical photons, whereas protons and atomic nu- 
clei can generate gamma rays via the decay of neutral 
pions produced in hadronic interactions with ambient 
material. A confirmation of the latter process associated 
with the gamma rays from SNRs is seen as a probe of 
nucleonic cosmic ray production sites in the Galaxy. 

The SNR G54. 1-1-0.3 has been predicted to emit 
gamma rays detectable with the current generation 
of TeV instruments after the discovery of a rotation- 
powered puls ar in its centre, a source of high energy 



electrons (e.g. Bcdnarek & Bartosik 

sar, PSR J1930+1852 (P^ 

down power oi E — 1 .2x10^^ erg /s and a spindown 



pOOSp . This pul- 
=136 ms), with a derived spin- 



age of rc~2900 years (Camilo et al 



2002), 



sars 



among 
known. The 



the youngest and most energetic pu 

PWN powered by this pulsar was first reco gnized as 

suc h from high resolution radio observations (Reich et 



al. 1985). The PWN exhibits a filled-center morphol- 
ogy (2' X 1.5'), a flat radio spectrum (a « 0.13) which 
is fairly constant over t he entire source extent, a nd a 
high level of polarization ( Velusamy fc Becker|1988 [). The 
discovery of an X-ray jet, consistent with the radio ex- 
tension, and the measurement of its non-thermal spec- 
trum confirmed the PWN nature of G54. 1-1-0.3 (Luet 



aOSOOT). More detailed X-ray observations by Lu et al 



] 2002 1 resolved the inner structure of the nebula as a cen 
tral point-like source with a surrounding ring and bipolar 
elongations, comparable to the Crab Nebula. The same 
authors suggested a distance of ~5 kpc, based on the X- 
ray absorption. The most recent estimate is 6.2^Q'g kpc, 
based on an apparent morphological associatio n between 
the PWN and a nearby molecular CO cloud (Leahy et 
al 



g 



2008) 



Similar to the Crab Nebula, G54.1-I-0.3 lacks a strong 
surrounding circumstcUar interaction between t he SNR 
forward sh ock and the progenitor star wind (Cheva- 
lier ||2005|), in which particles could be accelerated and 



produce an ad ditional non- t herma l emission contribu- 



tion. Recently, Leahy et ah] (2008) found a shell of 30' 
diameter surrounding G54.1-f-0.3, but the authors con- 
cluded that these are unlikely to be associate d based 
on the di fferent measured distances. After that, Koo et 
al. ( 2008[ ) identified a partial shell-like infrared structure 
surrounding much of the X-ray/radio core in G54.1-f 0.3, 
which they proposed to be a star-forming region triggered 
by late-phase winds from the SNR progenitor star. Based 

^^ TeVCat. an onlin e catalog for TeV Astronomy, 
|http://tevcat.uchicago.edu/| 



on hig h-resolution IR spectra, however, |Temim et al.| 
({20091 suggested that the shell is composed of freshly- 
formed ejecta dust that appears to be heated by the ex- 
panding PWN. 

G54.H-0.3 has not been previously detected at 7- 
ray energies. The source lies 1.4° from the EGRET 
source 3EG J1928-fl733 ( [Hartman et~al||1999[) . This is 
slightly outside the 99% confidence-level error box. The 
older (rc~8.2xl0^ years) and less-energetic {E = 
1.6 xlO^^ ergs/s) pulsar PSR J1928-)-1746 is proposed 
as a m ore likely coun terpart to the unidentified 3EG 
source (Cordes et al. 2006). This is strengthened with 
the much better localization of the gamma-ray source 
IFGL J1929.0+1741C, asso ciated with 3EG J1928-f 1733 
in the first Fermi catalogue (|Abdo et al.|2010[ ) . This other 
pulsar is 1.2° away from PSK J 1930-1-1852. In the VHE 
7-ray domain, the only published observation is an upper 
limit from the HEGRA collaboration at 20% of the Crab 

In 



Nebula flux above 600 GeV (Aharonian et al 



2002) 



this Letter we report the discovery by V ERl'l AS of 
VHE gamma-ray emission from the direction of the very 
young system G54.H-0.3/PSR J1930-hl852. The instru- 
ment and observations are described in § 2, while the 
analysis performed on these data is detailed in § 3. Re- 
sults are presented in § 4 and interpretation and conclu- 
sion in § 5. 

2. OBSERVATIONS 

VERITAS is an array of four imaging atmospheric 
Cherenkov telescopes located at the Fred Lawrence 
Whipple Observat ory in southern Arizona, USA ( Holder 



et al. 2006 2008). Each telescope has a mirror area ot 



110 m^ and is equipped with a 499-pixel camera of 3.5° 
diameter field-of- view. The system, completed in the fall 
of 2007, is run in a coincident mode requiring at least two 
of the four telescopes to trigger in each event. This design 
allows for the observations of astrophysical objects in the 
energy range from 100 GeV to above 30 TeV. VERITAS 
has an angular resolution (68% containment) of 0.1° di- 
ameter per event and a 5cr point source sensitivity of 1% 
of the steady Crab Nebula flux above 300 GeV in less 
than 50 hours of observation at a 20° zenith angle. 

The first pointed observations of the G54.1-F0.3/PSR 
J1930-I-1852 system with VERITAS took place during 
a survey of northern-hemisphere Galactic pulsars be- 
tween September 2007 and October 2008. The live time 
of this first data set, after data-quality selection based 
on weather and hardware status, amounts to 14.5 hr 
(some of which was taken in partial moonlight). Anal- 
ysis of these data revealed evidence of a VHE 7-ray sig- 
nal at the position of PSR J1930-I-1852 at the Aa level 
after accounting for all analysis trials. This evidence trig- 
gered further observations between April 2009 and June 
2009, resulting in 22.2 additional hours after quality se- 
lection. Of these, the most recent 15.2 hours of data were 
taken using only three telescopes, as telescope Tl was in 
the process of bein g relocated to optimize the instrument 
sensitivity (Perkins et al. 12009). The zenith angle range 
of observation of the combiiied data sets is 13° to 35°, 
and the total live time is 36.6 hr. 

3. ANALYSIS 

After calibration, the standard VERITAS event recon- 
struction scheme is applied to the data: tail-cuts image- 



TABLE 1 

Analysis results at the high Edot pulsar locations in the FOV of the PWN G54. 1+0.3 



PSR name 


Offset 

n 


On 

Events 


Off 

Events 


a^ 


Excess 
Events 


LiMa 
Significance 


Flux >1 TeV 
[10-i^cm-2s-i] 


J1930+1852 
J1928+1746 


0.5 
0.7 to 1.7 


231 
108 


720 
509 


0.18 
0.19 


101 
13 


7.0(7 
1.2 a 


5.3 ± 1.3 
<2.6 



Note. — The analysis used for the results presented in this table corresponds to an image integrated charge '-~~'225 photoclectrons and Q < 0.015 
dcg" 
^ Normalization factor for tlic different aeccptance of the background and source regions in the ring background model 

clean ing, second-m onient parameterization of these im- 
ages ( Hillas|[T985 ) and stereoscopic direction determina- 
tion based on ttie i nterse ction of image axes (see e.g. 



Krawczynski et al. (20061 for details). An initial selec- 
tion ot the events is performed to reject those in which 
the primary parameters, energy and direction, cannot 
be well-reconstructed. The events surviving this quality 
selection are those for which at least two telescope im- 
ages have more than 4 pixels included in the image with 
a total integrated charge in each image of more than 
80 photoclectrons, and with less than 20% of the total 
charge in the outermost ring of camera pixels. This se- 
lection removes images with insufficient Cherenkov light 
for accurate parameterization, or which are truncated by 
the camera boundary. In addition, events for which only 
the telescopes Tl and T4 have surviving images are ex- 
cluded, as their separation of only 35 m introduces an 
irreducible high background rate, du e to local muons 



a> 



deg rading the instrument sensitivity (Maier & Knapp 
12007^ 

TTTurther selection of events is performed to suppress 
the otherwise overwhelming hadronic cosmic-ray back- 
ground. This is achieved by comparing the shapes of 
the event images in each telescope with the expected 
shapes of simulated gamma-ray showers. These mean- 
reduced- scaled width (MRSW) and me an-reduced-scaled 



leng th (M RSL) cuts (see definition in Krawczynski et 
al. (2006)), as well as an additional cut on the ar- 
rival direction of the incoming gamma-ray with respect 
to the source position {9'^) reject more than 99.9% of 
the hadronic cosmic-ray background while keeping 45% 
of the gamma-rays. Different sets of these cuts, op- 
timized for different source properties, are applied to 
the G54. 1-1-0.3 data. Two a priori test positions are 
defined as the locations of the high spindown power 
pulsars: J1930-t-1852 (19''30'"30.13^ -|-18°52'14.1") and 
J1928-I-1746 (19'*28"42.48",-|-17°46'27"). In the search 
for a VHE 7-ray signal from a point source with a Crab- 
like spectrum the following cuts were employed: -1.2 < 
MSCW/L < 0.5 and 6"^ < 0.015 deg^, resuhing in an 
analysis threshold of 250 GeV. Since the emission region 
may be extended, as is the cas e of most of the galacti c 



PWN detected at VHE 7-rays ( [Aharonian et al.||2005[ ), 
a larger cut in the arrival direction 9"' < 0.055 deg^ is 
also considered. In addition, a rather tight cut on the 
integrated charge of ~225 photoclectrons is applied, to 
maximize the signal-to-noise ratio for a weak source with 
a hard spectrum. The latter cut results in an analysis 
threshold of 500 GeV and in superior angular resolu- 
tion. For the spectral analysis, this cut is 80 photoelec- 
trons providing maximal coverage in energy. 
For the background estimation, we employ two differ- 
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Fig. 1. — Smoothed excess map from G54. 1+0.3, as measured 
by VERITAS. Only events with image sizes above 225 photoclec- 
trons are used. The color scale indicates the number of excess 
gamma-ray events. The high E pulsars present in this field-of- 
view, J1930+1852 and J1928+1746, are marked with white plus 
signs. The contour of the 99% confidence level error box of the 
unidentified EGRET source 3EG J1928+1746 is overlaid in black. 
The gamma-ray sources in the first Fermi catalogue are marked 
with black crosses . 



ent methods, as described in Berge et al. (2007). The 
background estimate for each position in the two- 
dimensional sky map is taken from a ring of mean ra- 
dius 0.5° and an area six times that of the on-source re- 
gion. For the spectral analysis, the background is taken 
from 7 positions in the field-of-view with the same offset 
from the pointing direction as the source regions. This 
eliminates the need for corrections concerning the radial 
dependence of the background acceptance. The on-source 
and background region counts, together with a normal- 
ization factor for the different acceptance of these re- 
gions, are used to derive the statistical sign ificance of 
any ex cess following the likelihood method of |Li fc Ma | 
(fT983|. 



4. RESULTS 

Figure [Tl shows a sky map of the excess counts for 
the entire aata set in the region around G54.1-I-0.3/PSR 
J1930-I-1852 system, derived using the point source anal- 
ysis described above. A source of VHE gamma rays is 
clearly visible and is coincident with the position of the 
X-ray and radio PWN. An excess of 101 events is ob- 



served, corresponding to a pre-trials statistical signifi- 
cance of 7.0a at the location of PSR J1930-hl852, as 
shown in Table [l] This significance is the highest ob- 
tained at the pulsar position for the four different search 
criteria attempted, and this is taken into account in as- 
sessing the true chance probability of the observed signal 
of 6.8ct. As discussed above, the earlier data set showed 
evidence for a VHE 7-ray signal at 4.3ct. This evidence 
is confirmed by the presence of a 5.5cr detection at the 
same position in 2009 data set alone. This increase is 
as expected for a steady source with increased exposure 
time and a slightly reduced sensitivity of the array with 
only three telescopes. 

The unsmoothed map of excess events, binned in 
0.05° xO. 05°, is fit by a two-dimensional Gaussian in 
order to study the morphology of the observed emis- 
sion. The gamma-ray excess is well fit by the point- 
spread function of the instrument at the correspond- 
ing analysis threshold of 500 GeV (^0.11°), which re- 
sults in a detection compatible with a point-like source 
for VERITAS. The best fit of the Gaussian centroid is 
19''30"32" ± 25^+18°52'12" ± 20" (J2000) and hence 
we assign the name VER J1930-I-188. Given the system- 
atic uncertainty of 0.02° in this measurement, the cen- 
troid is consistent with the pulsar position (at a distance 
of 0.007°) and the 0.03° extent of the PWN G54.1-t-0.3 
around the pulsar. 

The total 7-ray excess above 250 GeV, i.e. that within 
a circle of 0.15° around the pulsar position, is 214±43 
events. Figure [2] shows the reconstructed gamma-ray 
spectrum from these events. The data are consistent with 
a power law in energy that extends from 250 GeV up to 4 
TeV with a photon index of T = 2.39 ± 0.23stat ± OM^ys 
and a differential flux at 1 TeV of {7 .5ztl.2stati^-5sys) x 
-[^Q-13 XeV~^ cm~^ s~^. The corresponding flux above 1 
TeV is 2.5% of the Crab Nebula flux. 

No other signiflcant TeV source is found in the maps, 
including at the second a priori test position, the loca- 
tion of PSR J1928-I-1746, for any of the searches per- 
formed. The signiflcance at this pulsar position is 1.2cr 
and the flu x upper limit above 1 TeV at the 99% confi- 
dence level ( Feldman fc Cousins [1998 1 assuming a power 
law source spectrum with P = 2.5 is F < 2.6x 10~^^ cm^^ 
s"""^ (about 1.4% of the flux of the Crab Nebula; see Ta- 
ble fl|. We have used the point source analysis for these 
calculations, as there is no known X-ray or radio PWN or 
SNR counterpart associated with this pulsar which can 
be used to constrain the size of the TeV emission. 

5. DISCUSSION AND CONCLUSION 

VERITAS has detected VHE gamma-ray emission 
from the direction of the SNR G54. 1-1-0.3, providing 
a new input to study the non-thermal particle pop- 
ulation in this object. This new TeV source, VER 
J1930-I-188, appears point-like within the instrument res- 
olution and does not show a notable offset from the pul- 
sar PSR J19304-1852, placing it among the more compact 
TeV sources found to be associated with an SNR. The 
emission is co-located with the X-ray and radio PWN in 
the SNR G54.1-I-0.3, and also with a molecular CO cloud 
in apparent morphological association with the PWN 
and a partial shell of IR ejecta around the PWN. 

Assuming the gamma-ray production in VER 
J1930-I-188 is related to the pulsar, the compactness of 
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Fig. 2. — Reconstructed VHE 7-ray spectrum of the new source 
VER J1930+188. The 99% confidence level are determined with 
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the TeV source and lack of significant offset between 
the pulsar and its X-ray PWN and the TeV source 
can be attributed to its relative youth. The spin-down 
age of PSR J1930-H1852 is 2900 years, which places 
it as the latest addition to an emergent class of TeV 
sources associated with very young pulsars: PSR 
B0531-h21/Crab Nebula ( |Weekes et al l [l989l), PSR 
J1747-2809/HESS J1747-281 i n G0.9-H0.1 (Camilo et al] 



1 |2009| [Aharonian et al.||2005[ ), PSR J1845-0258/HES 
7TM?-A;j9 in SiNR Kes 75, PSR J1833-1036/HESS 



J1833-105 in G21.5-0.9 ( Djannati-Atai et al 
PSR J1813-1749/H E5H 



20071 



and PSR J1813-1749/HESS J 1813-178 m 
0.02 (Gott heff k Halpern|[2009| . Those PWNe are good 
candidates tor being in the early stage of evolution, 
where the wind nebula is interacting with the freely- 
expanding supernova ejecta. This population stands 
in contrast to an apparently more common variety of 
energetic-pulsar-associated TeV sources, a group which 
contains, a s notable examples, PS R B1823-1313/IIESS 
J1 825-137 dAharonian et al.|2006a[ ), PSR B0833-45/Vela 
X ([Aharonian; et al.||2006b| andPSR J2229-F6114/VER 
J2227-h608 (lAcciari et 'all l2009[ ). These objects show 
extended TeV emission which is both larger than the 
X-ra y PWN and off set fro m the pulsar location (see 
e.g, de Jager & Djannati-Atai ([2008)) 



A 



common 



element in this class is older, middle-aged pulsars, with 
ages >10^ yr. Here the PWN may be displaced in some 
direction, possibly by an asymmetric SNR reverse shock 
generate d from an explosio n into an inhomogeneous 



medium (Blondin et al. 



2001). 



Table 2 lists the pulsar properties (spin-down power, 
age and distance) and the nebular energetics in several 
compact PWNe detected in TeV. The total power, Lg, 
radiated by G54.1-h0.3 in the energy band 0.2 - 10 TeV 
is 2.1 X 10^'' erg s~^ compared to 8.6 x 10^** erg s~^ 
for the Crab Nebula. This value of the gamma-ray lu- 



TABLE 2 

Pulsar and nebular properties of compact PWN detected at TeV 



Name 



E 
(10^'' erg/s) 



Tc 

kyr 



Distance 
(kpc) 



(10'="' erg/s) (10^'' erg/s) 



,34 



(% ) (% ) 



Crab Nebula 


460 


1.2 


2 


8.6 


1900 


0.02 


4 


220 


GO. 9+0.1 


43 


5.3 


8.5 


4.4 


36 


0.1 


0.8 


8 


G21. 5-0.9 


33 


4.85 


4.7 


0.72 


22 


0.02 


0.7 


30 


G54. 1+0.3 


12 


2.9 


6.2 


2.1 


2.9 


0.17 


0.2 


1.4 


Kes75 


8.1 


0.7 


7.5 


2.5 


24 


0.3 


2.9 


10 



Luminosity in tlic 0.2-10 TcV band 

Luminosity of the PWN in the 0.5-8 keV 

Efficiency to convert the spin-down energy 

Efficiency to convert the spin-down energy to X-ray synchrotron : rjx=^x/E 



band, taken from [ Kargalstev & Pavlov | l| 2008[ l 
/ to VHE gamma-ray emission : rjry=h~y /E 



minosity of G54.1+0.3 is the same order of magnitude 
what has been found for the other PWNe in the ta- 
ble. A more relevant quantity to compare is the efficiency 
of TeV luminosity production compared to the pulsar 
power which is rj^ ~0.17% for G54. 1-1-0.3. This is similar 
to that of GO.9+0.1 (« 0.1%) and 10 times larger than 
other young systems like the Crab Nebula (~ 0.02%) and 
HESS G21.5-0.9 (w 0.02%). However, ah these PWNe 
have significantly lower levels (ranging from 0.1 to 2.9%) 
of X-ray luminosity production efficiency than the Crab 
Nebula (« 4%). Despite being worst synchrotron emit- 
ters than the Crab Nebula, G54.H-0.3 together with the 
other recent compact TeV PWNe, are similar and even 
better TeV emitters than the Crab Nebula. 

Another diagnostic is the strength of the synchrotron 
emission relative to the gamma-ray emission, L2:/L^. The 
Crab Nebula has the largest ratio among all young sys- 
tems by far with L^/L^ of 220. That of G54.1-h0.3 is 
the lowest among the others that range between 1.4 and 
30. These numbers are consistent with 'Mattana et. al.l 



( 2009 1 who finds L^ >Lx after 5 kyr from pulsar birth, in- 
cluding a different list of objects and spanning over a 
larger range of Tc- The scattering of this ratio in two or- 
ders of magnitude, may be attributed to the unique envi- 
ronments in each of these objects, such as the local ambi- 
ent photon fields or the interactions with the surrounding 
mediums resulting in a n increase o f the n e bular magnetic 



field. As su ggested by Chevaher (2004); Djannati-Atai 



et al. (2007), these PWJNe may not be in energy equipar 



tition between particles and magnetic fields, as is the 
case for the Crab Nebula. Instead, these may be particle- 
dominated objects. However, proper modeling of each of 



these PWNe and its surroundings is needed to constrain 
the nebular magnetic fields. 

(|2008|) of the 



al. 



Given the association by Leahy et 
PWN with a molecular CO cloud, an alternative inter- 
pretation of the TeV emission from G54. 1-1-0.3 would be 
radiation from the decay of neutral pions produced in 
the hadronic interactions produced between the dense 
target and the SNR. To establish if this is a compati- 
ble scenario, deeper observations of this source would be 
needed to constrain the size and location of the acceler- 
ated particles in the TeV range. A detection by Fermi- 
LAT would add data to the spectral energy distribution 
of the system, giving an idea of the possible role played 
by protons by showing (or not) the typical bump in the 
spectrum exp ected from neutral pions (e.g. AharoniaiiJ 
Drury fcVolk ( 1994) ). However, based on the association 
with a high spin-down pulsar, a comparable size of the 
X-ray and radio PWNe given the instrument resolutions 
and comparison with similar detected objects, the PWN 
scenario best explains the data. Under the assumption 
that the gamma-ray production in VER J1930-I-188 is 
IC emission, the recently discovered IR shell of ejecta is 
an additional source of seed photons from the local dust 
up-scatted into the TeV energies and contributing in the 
observed emission. 
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